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ABSTRACT 


The  design  of  filters  and  amplifiers  with  minimum  element  distribu¬ 
tion,  minimum  number  of  elements,  minimum  sensitivity  to  parameter  change, 
etc.  is  important  today.  Since  the  synthesis  problem  usually  has  an 
infinite  number  of  solutions,  it  is  desirable  to  be  able  to  transform  a 
given  design  into  a  network  with  the  same  transfer  function  but  having 
fewer  elements  or  better  element  distribution,  etc.  The  techniques  of 
equivalent  network  theory  whidi  apply  a  congruence  transformation  to  a 
network  matrix  in  order  to  generate  equivalent  networks  are  used  to  do 
this.  An  extension  of  the  Cauer  formulation  is  discussed  which  applies 
the  transformation  directly  to  the  element  parameter  matrices,  thereby 
placing  the  elements  of  the  transformed  network  directly  in  evidence. 

The  transformation  is  then  used  to  transform  to  an  eqiri valent  network 
whose  elements  differ  from  those  in  the  original  network  by  only  an  in- 
crementail  amount.  In  the  limiting  case,  a  set  of  differential  equations 
for  the  elements  of  the  equivalent  network  result  with  arbitrary  inputs. 
The  method  of  steepest  descent  is  applied  to  choose  these  inputs  in  order 
to  force  the  network  to  converge  on  the  optimal  design. 


OONTINDOUSLY  EQUI7A11NT  NETWORKS  AND  THEIR  APFLICATION 

By 

J,D,  Schoeffler^ 


I,  Introduction 

Among  the  major  goals  of  network  theory  are  the  elimination  of 
transformers  from  networks  and  the  control  of  element-value  distribution, 
topology,  and  sensitivity  of  the  network  to  changes  in  element  values,^ 
Usually  these  goals  are  approached  by  demanding  that  a  synthesis  technique 
be  found  which  yields  a  network  with  the  desirable  qualities  listed  above. 
Hence  a  great  nvmber  of  different  synthesis  techniques  have  been  found,  no 
one  of  which  is  general  enough  to  attain  all  of  the  listed  goals.  For 
example,  the  Brune  and  Darlington  synthesis  methods  yield  networks  with  a 
minimum  munber  of  elements  but  have  mutual  coupling  among  the  elements. 

The  Bott-Duffin  method  requires  no  transformers  but  uses  an  extremely  large 
number  of  elements. 

A  different  design  philosophy  was  propounded  by  Cauer  in  1929  and 

3 

Howitt  in  1931,  They  found  that  a  whole  family  of  networks  could  be 
generated  from  a  given  network  by  means  of  matrix  transformations.  If 
these  transformations  were  properly  selected,  all  of  the  networks  so  found 
were  equivalent,  that  is,  had  the  same  driving  point  and/or  transfer 
impedance  as  the  original  network.  This  result  implies  that  the  great 
burden  of  achieving  all  of  the  listed  desirable  goals  can  be  removed  from 
the  synthesis  technique.  The  synthesis  method  need  merely  find  a  single 
suitable  network.  Then,  from  the  group  of  equivalent  networks  found  by 
matrix  transformation,  a  network  is  selected  which  has  no  mutual  coupling, 
and  does  have  good  element  distribution,  topology,  etc. 


•1 


I 


-2- 


Unfortunately,  the  matrix  transformation  method  never  realized  its 
potentiality  despite  much  effort.  *  *  The  problem  is  that  there  is  no 
unique  relation  between  network  topology,  elements,  and  parameter  matrices. 
Hence,  as  soon  as  the  original  parameter  matrix  is  transformed,  it  becomes 
very  difficult  to  recognize  the  new  network.  In  particular,  it  is  diffi¬ 
cult  to  see  in  advance  how  to  choose  the  transformation  matrix  so  that 
desirable  elements  result.  Only  the  simplest  applications  of  this  method 
have  thus  far  been  found. 

In  this  paper,  an  extension  of  the  Cauer  formulation  of  the  equivalent 
network  problem  is  discussed  which  applies  the  transformations  directly  to 
the  element  parameter  matrices,  thereby  placing  the  element  of  the  trans¬ 
formed  network  directly  in  evidence.  The  transformation  is  then  used  to 
transform  to  an  equivalent  network  whose  elements  differ  from  those  in  the 
original  network  by  only  an  incremental  amount.  In  the  limit  as  this 
Incremental  amount  approaches  zero,  a  set  of  first  order  linear  differential 
equations  for  the  elements  of  the  equivalent  network  result.  The  elements 
of  the  original  network  then  serve  as  the  initial  conditions  for  these 
equations .  The  differential  equations  have  some  arbitrary  input  functions 
such  that  for  any  choice  of  these  inputs,  the  solutions  to  the  equations 
yield  elements  which  are  functions  of  the  independent  variable  and  which 
for  any  value  of  the  independent  variable  constitute  an  equivalent  network. 
Thus  it  is  possible  to  transform  from  a  given  network  to  an  equivalent 
network  by  varying  the  elements  of  the  original  network  in  a  continuous 
manner.  The  advantage  of  this  approach  is  that  the  elements  of  the  trans¬ 
formation  matrices  become  very  easy  to  choose  in  contrast  to  the  situation 


where  the  transformation  is  not  continuous 
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The  problem  of  selecting  these  arbitrary  inputs  so  that  optimal 
networks  result  is  considered.  In  particular,  the  method  of  steepest 
descent  is  applied  to  design  networks  with  minimum  element  distribution 
and/or  minimum  ntiraber  of  elements. 

II,  Equivalent  Network  Theory 

7  8  9  10 

Recently  it  has  been  shown  *  *  *  that  the  Cauer  formulation  of 
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the  equivalent  theory  problem  could  be  extended  in  such  a  way  that  direct 
control  of  the  elements  of  the  equivalent  network  resulted.  To  be  more 
explicit,  the  usual  Cauer  formulation  generates  equivalent  networks  by 
applying  conginience  transformations  to  the  terminal  admittance  matrix  of  a 
network.  This  technique  was  extended  so  that  transformations  can  be 
applied  directly  to  the  branch  admittance  matrix  while  at  the  same  time 
preserving  a  driving  point  or  transfer  ratio  or  some  combination  of  driving 
point  and  transfer  ratios.  These  results  are  used  in  this  paper  and  are 
summarized  here. 

let  the  equilibrium  equations  of  a  network  with  n  independent  mode  pairs 
on  the  admittance  basis  be  given  by 


I  -  E  2.1 

where  I  is  the  nxl  column  matrix  of  current  variables  and  E  is  the  nxl  column 
matrix  of  voltage  variables  and  is  the  nxn  terminal  admittance  matrix. 

To  transfom  to  an  equivalent  network,  Cauer  applies  the  transformation: 

E  “  A.  E' 

4 


2.2 
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\rtiere  is  nonsingular  and  superscript  t  means  transpose.  This  results 
in  the  new  set  of  equilibrium  equations 


I'  =  Y'  E' 
4 


2.3 


where  the  new  terminal  admittance  matrix  is  related  to  the  given  matrix 
by  a  congruence  transformation: 


Yj,  =  Y.  A, 
4  4  4  4 


2.4 


Equivalence  is  assured  if  A^  is  nonsingular  and  certain  elements  of  A^  are 
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properly  chosen.  That  is,  the  new  and  old  networks  are  assured  of  being 
reali2able  and  of  having  any  combination  of  driving  point  and/or  transfer 
impedances  the  same. 

By  using  the  generalized  equilibrium  equations  of  Guillemin,  it  was 
shown  that  a  related  transformation  could  be  applied  directly  to  the  bxb 
branch  admittance  matrix  Y  where  Y,  is  defined  by 


J  -  Y,  V  2.6 

D 

and  J  and  V  are  the  (bxl)  column  matrices  of  branch  voltages  and  currents. 
Making  the  transformations 


V  -TV 


J'  -  T''  J 


yields  a  new  branch  admittance  matrix  Y^  given  by 


Y '  -  T  Y^  T 
b  b 


2.6 


2.7 
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Since  the  elements  of  and  Y^  are  the  elements  of  the  original  and 
transformed  networks,  it  is  clear  that  this  approach  to  equivalence  gives 
the  elements  of  the  new  network  explicitly,  thereby  allowing  control  over 
them.  Note  also  that  if  the  original  network  is  passive,  the  branch 
admittance  matrix  is  positive  definite  or  semi -definite  and  that  the  con¬ 
gruence  transformation  preserves  this  characteristic,  If  the  original  net¬ 
work  has  no  mutual  coupling,  Y^  is  a  diagonal  matrix  with  all  positive 
entries.  If  T  is  chosen  so  that  Y^  is  also  diagonal,  the  congruence  trans¬ 
formation  insures  that  all  of  the  elements  of  Y^  are  also  positive.  Thus 
no  negative  elements  can  appear  due  to  the  application  of  the  transformation. 

5 

This  is  the  problem  in  the  usual  Cauer  formulation  as  discussed  by  Howitt. 

The  most  general  branch  is  a  parallel  connection  of  R,  L,  and  C  and 
hence  the  branch  admittance  matrix  of  the  network  may  be  written  as 

Y.  *=  G  +  sC  +  i  2.8 

b  s 

where  G  is  the  branch  conductance  matrix,  C  is  the  branch  capacitance  matrix, 
and  L  ^  is  the  branch  inverse  inductance  matrix.  It  has  been  shown  that 
equivalent  networks  result  if  congruence  transformations  are  applied  to  each 
of  the  element  matrices  as  follows: 


G '  ■=  T  “  G  T^ 


C  =  Tj  C  Tg 


(L'^) '  =  Tj  (L'^)  Tj 


2.9 


and  if  each  of  the  transformations  satisfies  the 
equation 


Ti 


A. 


following  constraint 


2.10 
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where  a  is  the  usual  nxb  cut  set  matrix  of  the  network  and  is  a  non¬ 

singular  nxn  matrix  with  certain  rows  properly  chosen  to  keep  desired 
transfer  ratios  invariant.  Notice  that  different  transformations  may  be 
applied  to  each  of  the  element  matrices  but  all  must  satisfy  Eq.  2.10 
with  the  same  non singular  A^. 

If  the  original  network  was  realizable,  the  transformed  network  is 
also  realizable  since  the  congruence  transformation  leaves  positive 
definite  and  semidefinite  matrices  at  least  positive  semidef inite .  In¬ 
spection  of  Eq.  2,7  'shows  that  the  new  network  defined  by  has  no 
mutual  coupling  provided  that  TtY^T  is  diagonal  or  that,  in  other  words, 
each  of  the  is  orthogonal  with  respect  to  the  corresponding  parameter 
matrix.  The  placing  in  evidence  of  the  elements  of  the  new  and  old  net¬ 
works  in  this  approach  to  equivalence  is  very  advantageous.  The  problem 
of  choosing  the  arbitrary  elements  of  the  transformation  in  order  to 
arrive  at  optimal  networks  is  the  subject  of  this  paper, 

III.  Continuously  Equivalent  Networks 

The  transformation  from  one  network  to  another  equivalent  network 

can  be  visualized  by  considering  the  vector  space  interpretation  of  a 

network.  Let  a  network  have  b  elements.  Then  the  branch  admittance 

2 

matidx  Y^  is  square  with  b  entries.  For  a  reciprocal  network,  Y^  is 

z 

symmetrical  and  hence  there  are  only  (b  +  b)/2  independent  quantities  in 
the  matrix.  Call  this  number  Q.  In  a  Q-dimensional  vector  space,  each 
coordinate  of  vrtiich  is  the  value  of  one  entry  of  Y^,  each  point  represents 
a  network.  The  set  of  all  networks  equivalent  to  a  given  network  corres¬ 
ponds  to  a  set  of  points  on  a  surface  in  this  Q-dimensional  soace.  If  an 
equivalent  network  has  no  mut\ial  coupling,  only  b  of  the  Q  entries  are 
nonzero  and  the  set  of  all  equivalent  networks  without  mutual  coupling 


I 


-7- 


corre spends  to  a  set  of  points  on  a  subsurface  of  the  surface  of  equiva¬ 
lent  networks.  We  are  interested  in  transforming  from  the  original 
network  (a  point  in  Q-space)  to  a  more  desirable  network,  that  is,  to 
some  other  point  in  Q-space  which  lies  on  the  surface  of  equivalent  net¬ 
works. 

For  simple  networks,  it  is  possible  to  choose  transformation  matrices 
which  satisfy  the  equivalence  restrictions  and  which  result  in  networks 
which  are  more  desirable  in  some  sense.  But  for  a  network  of  even  moderate 
complexity,  it  is  not  easy  to  choose  the  free  parameter  of  T.  With  the 
vector  space  internretation  in  mind,  we  can  imagine  the  transformation  not 
as  a  single  "jump"  from  one  point  on  the  surface  to  another,  but  rather 
as  a  succession  of  transformations  T^,  Tg,  ....  each  of  which  transforms 
the  original  point  to  another  point  on  the  surface  a  differential  distance 
away.  In  the  limit,  the  transformation  T  becomes  continuously  varying  and 
the  point  on  the  surface  of  equivalent  networks  moves,  following  a  contin¬ 
uous  line  from  the  original  network  to  the  final  network. 

To  make  this  approach  more  concrete,  consider  an  independent  variable 
X.  Then  if  the  transformation  T  is  considered  as  a  function  of  x,  T(x) 
the  kth  element  of  the  equivalent  network  e^^  is  also  3  function  of 
X,  ej^(x).  That  is,  for  any  value  of  x,  the  set  e^(x),  eg(x),  ....,e^(x) 
is  equivalent  to  the  original  network.  Using  this  approach,  the  selection 
of  the  transformation  matrix  T  is  greatly  facilitated. 

To  formulate  the  problem  on  a  continuous  basis,  consider  the  trans¬ 
formation  matrix  T  which  transforms  I,  into  I.  ♦  AY.  whei-e  AY.,  is  incre- 

D  b  b  b 

mental  in  size: 


T  »  U  +  AT 


3.1 


-8- 


where  U  is  the  bxb  unit  matrix.  The  application  of  T  to  results  in 


given  by 


y,  .  yUy  „ 

D  D 


y^  +  Y.  AT  +  AtV  +  AT^y.  AT 
b  b  b  b 


3.2 


Define 


T  *  T^Ax 


3.3 


where  Ax  is  a  scalar.  Then  T  and  Y*  become 

b 


T  -U+AT*U+  TyAx 


Y^  =  y^  +  (Y^T^  +  TyY^)Ax  +  O(Ax^) 


3.4 


Writing  T  -  T(x),  Y^  »  >  and  Y^  =  Y^(x  Ax),  Eq.  3.4  becomes 


b  "b' 


Y.  (x  +  Ax)  -  Y.  (x) 


f  0(Ax  ) 

_ X»__  ■  y  T  +  T  Y  +  5.5 

Ax  b  1  1  b  Ax 


Passing  to  the  limit  yields  the  matrix  differential  equation 


-3-2  -  Y.T-  +  t2^  Y, 
dx  b  1  lb 


3.6 


Let  Y^  =  Ze^.J  . 


Since  Y^  is  symmetrical  Eq.  3.6  is  symmetrical  and  hence  correspond  to 

O 

only  (b  ♦  b)/2  equations  each  of  which  is  of  the  form 


de.  ,  t  i“ij2,  ....,b 

-  uj  p,  T  +  T;  Y  3^u.  3.7 

to  i  i-b  1  1  b-i  3  j  ^ 
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til 

lAiere  is  the  k  unit  vector  in  b  space. 

The  matrix  is  a  function  of  x,  T^(x),  and  must  be  selected  so 
that  equivalent  networks  result.  To  derive  the  equivalence  restrictions 
on  T^,  start  with  the  general  restriction 

To^  -  a*  A.  3.8 

4 


Writing 


T  -  U  +  T^Ax 


3.9 


and 


A.  *  U  +  BAx 
4 


3.10 


we  have 


+  T^a^Ax  ■  a^BAx 


3.11 


or 


3.12 


Thus  the  restrictions  on  T^  are  very  similar  to  those  on  T  except  that 
B  is  of  the  form  (considering  the  preservation  of  the  driving  point  and 
transfer  impedances  at  terminal  pairs  1  and  2) : 


B  - 


0 

0 

X 

X 


0 

0 

X 

X 


0  0 

0  0 

X  X 

X  X 


3.13 


idiere  the  x's  indicate  arbitrary  numbers.  That  is,  B  has  rows  of  zeros 

corresponding  to  the  terminal  pairs  in  question.  Writing  out  Eq.  3.12 

2  2 

yields  nb  linear  equations  relating  the  b  entries  of  T^  and  the  b  -  kb 
entries  of  B  (k  is  the  number  of  terminal  pairs  at  which  equivalence  is 
retained).  Thus  fewer  than  b  of  the  elements  of  T  are  Independent. 
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The  technique  is  best  illustrated  by  an  example. 

Consider  the  network  shown  in  Fig.  1.  We  can  transform  this  network 
by  considering  it  to  be  two  parallel  networks  each  with  the  graph  and 
tree  shown  in  Fig.  2.  For  this  grajii,  tree,  and  branch  numbering,  the 


standard  cut  set  matrix  is 


-  [.; :  3 


If  the  branch  admittance  matrices  are 


-gi  0  0 

Gr  ■  0  g^  0 


‘ll  ^2  ^5 


^  *  ^2  ^22  ^25 

^3  ^25  ^33 


y,  -  G  +  i 
b  s 


According  to  Section  2,  these  may  be  transformed  to  find  equivalent 
networks  by  congruence  transformations  and 


G'  -  t;  G  T 
G  g 


d'h'  ■  t‘  l-\ 


provided  that 


•11 


t  t 
T^a  *  a 


3,19 


Trtiere  A^  is  an  nxn  nonsingular  matrix.  Assuming  only  the  Input  Impedance 


is  to  be  preserredy  A^  must  have  the  form 


1  0 


®1  ^2 


3,20 


with  ag  nonzero. 

In  the  continuous  case  we  have 


i  ■ 


and 


.  T  T  *  T  -1 

'IT  ^  ^2  *  ^2  ^ 


3.21 


where 


-  o*B 


TgO^  ■  a^B 


and 


B  -  A .  -  U  - 
4 


3.22 


In  order  for  A^  to  be  nonsingular,  bg  must  not  equal  -1, 
Let 


^5  ^6 

"8  "^9^ 


3.25 
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and 


n 

t' 

y 


4 

H 

H 


3.24 


Then  Eq.  3.22  give  constraint  relations  among  the  t.  and  aitong  the  t’  . 

^  J  ^  J 

Specifically,  each  matrix  equation  yields  6  homogeneous  equations  in  11 
parameters.  Hence  only  five  of  each  set  are  independent.  Solving  for 
the  remaining  six  in  terms  of  these  five  yields  the  constraints: 


-ti  - 

"l-^2 

-  t. 


-t?  .  b^ 

S  ^2 


-t .  -  b^ 

“^7 

H"^2 


3.25 


The  differential  equations  for  the  G  matrix  are  subject  to  a  further 
constraint.  That  is,  all  of  diagonal  elements  of  G  must  be  zero  since 
it  corresponds  to  a  resistor  network.  The  same  need  not  be  true  for  the 
inductance  matrix  because  mutual  coupling  is  allowed  in  this  case.  Hence 
Eqs,  3.21  becomes 


d 


or 


h 

0 

LO 


0 

0 


GT^  ♦  T£G 


dgi 

dT  “  ^1^1 


3.26 


IT  " 
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IT  "  ^5^9 

and 

0  -  ♦  g^tg 

0  «  g^tg  4  gjt^ 

0  “  ggtg  4  gjtg  3.27 

Confining  Eq.  3,27  with  Eq.  3,25  yields 

dgi 

IF  “  ^1^1 

dgj 

IF  “  ■  ^2*4 
IF  *  *  V 

0  -  -  ♦  bj)  »  gjt, 

0  “  *  ^3^7 

°  “  ®2^4  ®3^'^1  ■ 

The  parameters  t^,  t^,  and  t^  may  be  eliminated  from  the  differential 
equations  using  the  last  three  algebraic  equations.  Solving  yields 

ggitj  -  gj(bj  -  bj)  -  gjb^ 

Sgjt,  -  g^Cbg  7  b^)  -  gjbj^ 


Jgjt,  ■  gj(bj  ♦  bj)  ♦  gjb^ 


3.29 
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Hence  the  differential  equations  for  the  resistors  become 

dgi 

•sr  ■  -  *>!>  - 

^  -  -  gjCbj  ♦  bj)  .  gjb^ 

^  ’  g5(bl  *  2bj)  3.30 

Notice  that  the  algebraic  equations  in  Eqs,  3.28  may  not  be  solvable 
for  b^  and  bg  in  case  any  element  is  zero.  For  example,  if  ■  0  for 
any  value  of  x,  we  find 


0  =  ggt^ 

0  =  ggt.^ 

0  -  ggt^  -  g3(b,  -  V  3.31 

or 

t,  •  »!  •  »!  ■  0 

Hence  in  this  case,  b^  is  no  longer  independent.  Several  degenerate 
cases  of  Interest  arise.  These  cases  and  the  additional  constraints 
they  imply  are  sximmarized  in  Table  1,  Inspection  of  Table  1  shows  that 
as  an  element  becomes  zero,  it  must  stay  zero  and  the  additional  con¬ 
straints  can  be  derived  by  setting  that  element  equal  to  zero  in  Eqs. 

3,30,  This  follows  because  the  equivalent  networks  are  derived  from 
congruence  transformations  and  an  element  becoming  zero  causes  a  decrease 
in  rank  of  which  cannot  then  be  increased  by  a  congruence  transformation. 
The  differential  equations  for  the  elements  of  the  inductance  matrix 
are  not  so  restricted  since  off-diagonal  elements  need  not  be  zero.  Hence 
there  are  6  differential  equations.  For  the  purposes  of  illustration, 
suppose  that  no  mutual  coupling  is  allowed.  Then  only  three  differential 
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equations  result  and  they  are  Identical  in  form  to  those  for  the  resistors. 
The  resulting  set  of  six  differential  equations  for  the  network  is: 


ot  m 

®l 

glCbj  - 

bi)  - 

■ 

or  f  Mm 
®2 

■  * 

b^) 

'  esb 

ff  1  at 

®3 

bi)  ^ 

'  *  b’ 

h.  *  * 

nil 

hjiCbj 

"  V 

■  bsb 

hi  m  , 

"22 

"  V 

■bsb 

h’  “ 

33 

^  V 

*  .  b^) 

3.32 

This  set  is  valid  only  as  long  as  no  element  is  zero.  As  an  element 
passes  through  zero,  additional  constraints  as  in  Table  1  must  be  applied. 
The  inputs  b^(x)  and  bg(x)  are  arbitrary  functions  except  that  b^  must 
never  equal  -  1  so  that  remains  nonsingular.  Continuing  the  example, 
suppose  the  elements  of  the  given  network  are  those  in  Fig,  3  and  suppose 
we  arbitrarily  select 


bj^Cx)  -  1 

bg(x)  •  0  3.33 

Solution  of  Eqs.  3.32  with  inputs  defined  by  Eqs.  3,33  is  easily  accom¬ 
plished  by  hand.  The  resulting  elements  as  functions  of  x  are  shown  in 
Fig.  4a,  The  element  values  are  also  plotted  in  Fig.  5  for  x  between  0 
and  1/3.  It  is  easily  verified  that  for  any  x,  the  original  and  trans¬ 
formed  networks  have  the  same  driving  point  impedance.  At  x  =  1/3,  h^^ 
becomes  zero.  Thus  inductance  number  1  becomes  infinite  and  the  con¬ 
straint  b^  *  0  must  now  be  applied.  To  continue  transforming  the  network, 
suppose  we  now  choose 
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b^(x)  =  0 

bgCx)  “  1  3.34 

The  resulting  elements  for  x  >  1/3  are  showi  in  Fig.  4b,  and  are  also 
plotted  in  Fig.  5.  The  inputs  selected  here  were  completely  arbitrary 
and  merely  included  for  Illustrative  purposes.  The  selection  of  the 
arbitrary  inputs  in  order  to  arrive  at  optimal  networks  is  discussed  in 
the  next  section. 

IV.  Optimal  Network  Design 

As  mnetioned  in  the  introduction,  optimal  network  design  -  that  is, 
the  design  of  networks  with  prescribed  driving  point  and/or  transfer 
ratios  plus  minimum  element  distribution  and  minimum  number  of  elements 
etc,  -  is  difficult  using  present  synthesis  techniques.  The  following 
approach  is  proposed.  After  a  network  with  a  desired  transfer  function 
has  been  designed  by  network  synthesis  techniques,  the  network  is  trans¬ 
formed  continuously  to  a  more  optimal  network  by  proper  choice  of  the 
arbitrary  inputs.  For  example,  if  the  element  distribution  is  to  be 
minimized,  the  inputs  can  be  selected  so  that  as  x  increases,  the  element 
distribution  decreases.  Such  a  procedure  requires  an  explicit  error 
criterion,  that  is,  an  explicit  measure  of  the  deviation  of  the  given 
network  from  the  desired  optimum.  Since  solution  of  the  differential 
equations  for  the  elements  is  most  efficiently  handled  on  a  digital  corn- 
outer,  there  is  in  general  no  particular  advantage  to  choosing  very  simple 
linear  criteria.  Two  examples  of  error  criteria  are  the  following.  If 
element  distribution  is  to  be  minimized  a  convenient  measure  of  the 
deviation  from  the  optimum  is  a  quadratic  form; 

^  ~  2  (g  -  g)^ 


4.1 
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>ihere  g  Is  the  average  of  the  elements: 

g  -  i  Z  4.2 

k*l 

Thus  9^  is  a  mean  square  measure  of  the  deviation  of  the  elements  from 
their  mean  and  if  ^  is  minimized,  so  is  the  element  distribution.  Since 
the  g-  are  functions  of  x,  (jl  is  also  and  hence  if  we  calculate  dj^/dx,  it 

iC 

is  a  function  of  the  elements  and  their  derivatives  and  consequently  a 
function  of  the  arbitrary  inputs.  These  inputs  can  then  be  selected  to 
reduce  9^(x)  as  quickly  as  possible.  This  is  called  the  "method  of  steepest 
descent"  and  is  an  example  of  a  minimization  technique  that  can  easily  be 
applied."  The  similarity  to  the  optimization  problems  encountered  in 
the  field  of  automatic  control  should  be  noted." 

Another  criterion  of  Interest  is  the  number  of  elements.  For  example, 
in  an  RL  network  as  in  the  example,  it  may  be  desirable  to  minimize  the 
number  of  coupled  inductors,  etc.  A  possible  measure  of  the  deviation 
from  the  optimum  for  the  former  case  is  a  weighted  sum  of  the  inductance 
values.  Since  all  the  inductances  are  positive,  the  weightings  can  be 
chosen  to  emphasize  any  particular  element  or  group  of  elements.  For  the 
example  problem  we  choose 

where  aj^  are  arbitrary  weighting  factors.  The  derivative  of  ^^(x)  with 
respect  to  x  becomes  after  substituting  the  expressions  for  the  derivatives 
of  the  elements: 


4.4 
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dearly  to  minimize  d(2^/dx  (to  maximize  rate  of  descent)  choose 

pci'll  (a,  -  >1  -  aj)  a  .  ajl]<  0 

and 

Choosing  the  inputs  according  to  Eq.  4.5  and  solving  the  differential 
equations  is  best  handled  on  a  machine.  In  this  case,  the  differential 
equations  are  replaced  by  finite  difference  approximations,  requiring 
that  the  inputs  be  limited  in  magnitude.  According  to  Eq.  4.4,  the  arbi¬ 
trary  inputs  should  always  be  at  their  limits  with  their  signs  determined 
from  Eq,  4.5.  The  solution  of  the  set  of  equations  with  these  inputs 
is  straightforward  and  rapid  on  a  computer  and  poses  no  problems.  Notice 
that  the  weighting  factors  are  arbitrary  and  free  to  be  changed  during 
the  course  of  a  solution  in  order  to  weight  certain  elements  more  heavily 
at  the  discretion  of  the  designer  or  to  facilitate  rapid  convergence  to 
the  optimal  network.  In  general,  any  error  criterion  may  be  used  that 
can  be  formulated  quantitatively  and  the  arbitrary  inputs  then  chosen  to 
transform  the  given  network  into  a  more  desirable  one.  The  advantage  of 
the  continuous  transformation  over  the  discrete  approach  is  that  the 
choice  of  the  arbitrary  parameters  becomes  very  simple.  These  techniques 
have  also  been  applied  to  active  circuits  but  will  be  discussed  elsewhere. 
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V,  Suiamary 

An  extension  of  the  Cauer  formulation  of  the  equivalent  network 
theory  problem  is  discussed  lAiich  applies  the  transformations  directly 
to  the  element  parameter  matrices,  thereby  placing  the  elements  of  the 
transformed  network  directly  in  evidence.  The  method  of  steepest 
descent  is  applied  to  choose  the  arbitrary  parameters  of  the  transforma¬ 
tion  in  order  to  derive  optimal  networks,  that  is,  networks  with  mininnim 
number  of  elements  and/or  minimum  element  distribution. 


I 
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Table  I 

Summary  of  Differential  Equations  for  the  Resistors 


Case  Additional  Constraints  Differential  Equations 


No  Element 
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dg-i 
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1  ^ 
dx 

-  V 

■  ^3^ 

•  Sl(»2 

♦  b^) 

dx 

dgp 

dx 
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-bi) 

*  g3(2V^) 

gl  “  0 
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